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ABSTRACT 






The XR-3 Loads and Motions Program, a digital computer 
simulation of the six degrees-of- freedom equations of motion 
for the XR-3 captive air huhble test craft, is subjected to a 
wide range of regular and irregular wave excitations. The 
operating characteristics of the simulated craft and the 
frequency response functions for the heave, pitch, and roll 
responses in ahead, abeam, and following seas are obtained 
for regular wave excitations. The ahead seas frequency response 
functions are verified with irregular wave excitation. 
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I. INTRODUCTION 



A . BACKGROUND 

Since 1973 work has been ongoing at the Naval Postgraduate 
School on a digital computer model that simulates the forces 
generated in and the motion experienced by the school's two 
ton XR-3 captive air bubble test craft. The origin of this 
model was the "SES Loads and Motions Program" developed for 
the U.S. Navy by Oceanics Inc. for the Bell 100 ton (100-B) 
surface effect ship /~Ref. 1 _ 7 . Since there were major design 
differences between the 100-B and the XR-3, Leo and Boncal 
/~Ref. 2 _7 converted the 100-B model to an initial representa- 
tion of the XR-3 test craft. 

In 197^ Finley /~Ref. 3_7 refined the bow and stern seal 
subroutines of the XR-3 model and developed new lift fan maps 
in order to improve the model's performance by providing more 
realistic air flow rates. At the same time Forbes /~Ref. 4J7 
measured various variables of the actual XR-3 test craft while 
in operation and compared them with the computed values from 
the XR-3 model. Forbes also introduced several modifications 
to various subroutines of the model. By the end of 197^- there 
were three versions of a program to simulate the XR-3. 

During 1975 Menzel /~Ref. 5_7 undertook a comparison of 
Finley's and Forbes's versions by looking at their pitch and 
roll characteristics. Based on the results of this study, 
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Menzel introduced program modifications to improve the com- 
puted pitch and roll behavior of the simulated XR-3. Menzel' s 
version is considered to be the best to date but there are 
still portions of the model that require study. 

Before additional modifications are made, however, it was 
deemed desirable to provide a baseline with which further 
versions could be compared. In addition it was noted that 
the majority of the past work had been performed while oper- 
ating the model in calm water conditions. Although the model 
includes a "sea-state generator," relatively few runs had been 
made simulating the XR-3 craft in the presence of waves. Two 
reasons account for this. Past work has been primarily con- 
cerned with developing the force interactions within the 
model where the response of the model to wave forcing actions 
was not desired. Second, a large amount of CPU time was 
necessary to simulate a few seconds of simulated run time. 

The response of the craft, however, to various wave inputs 
provides a commonly used baseline for future comparisons. 

B. OBJECTIVE 

The objective of this study was to investigate the dynamic 
and operational characteristics of Menzel' s model when sub- 
ject to many differing operating conditions. In particular 
it was desired to obtain the frequency response of the model, 
i.e., the response of such characteristics as heave, pitch, 
and roll as functions of wave frequency and amplitude. In 
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support of this primary objective it was desired to review 
the characteristics of ocean waves and the interaction of 
waves and craft in the XR-3 digital model. 
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II. FREQUENCY RESPONSE FUNCTIONS 

The fundamental principle of linear systems is that in 
the frequency domain the response of a system, Y(w), is equal 
to the excitation, X(w) , multiplied by the system's transfer 
function, H(w) : 

Y(w) = X(w)H(w) . 

Also |Y(w)| =|X(w)| |h(w)| 

/Y(w) = /X(w) + /H(w) 

where | | denotes the magnitude and /_ the phase and 
S y (w) = S x (w)| H(w )| 2 

where S (w) and S (w) represent either the power spectral 

x y 

density or energy spectral density of the excitation and 
response respectively. 

This concept, ingrained in the electrical engineer, is 
also used by naval architects and marine engineers to describe 
the motion of a displacement type vessel in a random sea and 
was first proposed for such use by St. Denis and Pierson 

Ref. 6_J7. The marine engineer refers to the transfer func- 
tion as the "frequency response function" and | H(w)| as 
the "Response Amplitude Operator," (RAO) . The excitation is 
the ocean, wave and the excitation frequency is given in terms 
of the encounter frequency (see Appendix A). Using this 
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concept, the vessel can be represented as a single input, 
multiple output system as shown in Figure 1. The frequency 
response functions thus describe the vessel in terms of sev- 
eral linear input-output relationships. There can be a sep- 
arate frequency response function for each degree of freedom. 
They can be obtained either theoretically or experimentally 
by such techniques as towing a model in regular waves in a 
towing tank and measuring the response of interest. 




H^w) 

H 2 (w) 



R^t) 

R 2 (t) 



n(t ) 



> Responses 



wave 

displacement 




R ± (t) 



H n (w) 




Frequency 

Response 

Functions 



Craft Input-Output Relationship 



Figure 1 
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The assumption that the excitation-response relationship 
is linear is not true for the SES craft since the hydrodynamic 
and aerostatic characteristics of the SES are inherently non- 
linear. However a theoretical frequency response model for 
the XR-3 was developed by the Aerojet-General Corp. /“Ref. 7_7 
using a four degree-of-freedom representation of the craft 
(responses in surge, heave, pitch, and roll). It was assumed 
that the response motions were small enough so that the devia- 
tion from the linear response was negligible for the motions 
in the vertical plane. 

The frequency response characteristics obtained from this 
model are reproduced in Appendix B and discussed further in 
Chapter IV. The predictions were compared to the results from 
i scale model towing tests and, in general, the experimental 
and theoretical frequency response information was found to 
show good agreement. 

The purpose of this study was to obtain the frequency 
response functions for the XR-3 Loads and Motions program, a 
time-domain model that simulates the craft in six degrees-of- 
freedom and provides time domain outputs of the craft's motions 
and forces for the on-bubble mode of operation. Since the 
model provides time-domain outputs, determination of the 
response characteristics was basically an experimental process 
consisting of running the model in a sinusoidal wave many 
times for different wave frequencies, amplitudes, and craft 
speeds. The desired responses were observed, measured, and 
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plotted as functions of the encounter frequency. The process 
and results are discussed in detail in Chapter IV. 

However, because the program provides for all six degrees- 
of-freedom, the inherent nonlinearities of the SES craft are 
incorporated into the model, particularly in the coupling of 
the various degrees of freedom. The generality of the fre- 
quency response obtained is therefore dependant on the signi- 
ficance of the nonlinearities. 
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III. WAVE SIMULATION IN THE MODEL 



The XR-3 Loads and Motions program has the capability of 
simulating the craft's motion in regular and irregular seas. 
The purpose of this section is to review briefly the charac- 
teristics of ocean waves and their modeling in the computer 
program. 

A. WAVE CHARACTERISTICS 

The outstanding characteristic of the open ocean is its 
irregularity which can be described by statistical methods 
on the assumption that a large number of regular (sinusoidal) 
'waves having different wavelengths, frequencies, directions, 
phase, and amplitudes are superimposed to form the randomly 
varying sea. Consequently a brief look at the properties of 
regular waves is necessary. 

It is common to characterize wave motion by the direction 
of the field displacement; waves are either transverse or 
longitudinal. In the transverse wave the field intensity 
oscillates in a direction perpendicular to the direction of 
wave travel, as in the electromagnetic wave. In the longi- 
tudinal wave the particles of the medium move in a direction 
parallel to the motion of the wave itself as in the sound 
wave. However ocean waves can be considered a combination 
of the motion experienced in both the transverse and longi- 
tudinal waves since the path of an individual water particle 
is circular as a wave moves by. 
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The properties of regular waves are well defined by 
gravity- free surface wave theories. It is assumed that the 
wave crests are straight, infinitely long, parallel and 
equally spaced and that wave heights are constant. The wave 
advances in a direction perpendicular to the line of the 
crests at a uniform velocity, C, often referred to as "cel- 
erity" to emphasize that it is the wave form rather than the 
water particles that advances. 

For convenience the regular wave is normally modeled as 
a sinusoid with the following properties: 



Wave Velocity: 


C = 


( = 
K 2 7t ’ 


g T - 

2 -n w 


(III-D 


Frequency : 


w = 


= 


2 tt C 


(III-2) 


Wavelength : 


K = 


2 -7T g 
2 




(III-3) 



w 



Wave surface profile: 

A(x, t) = A n cos - wt + £. ) (III-4) 

O A 

where w = angular frequency (radians/sec) 

= wavelength 

T = wave period 

g = gravitational constant 

x = distance along the direction of 
propagation from a defined origin 

£. = phase angle at the origin 

A Q = maximum wave amplitude 
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The slope of the wave surface is obtained by differentiating 
(III-4) : 



= ~ 2 r tA ° sinAf^-wt+M 

x A A 

The maximum wave slope is then: 

2^rA 360 A 

v = — radians = — degrees (III-5) 

A A 

An important property of the regular wave is the energy 
associated with it consisting of kinetic energy associated 
with the orbital motion of the water particles and potential 
energy resulting from the change in water levels in wave 
hallows and crests. The average energy per unit area of 
surface for a regular wave is /~Ref. 8_7: 

2 

E = | k A q ( III-6) 

where k is a water density constant. 

The irregular elevation or amplitude of the sea surface is 
considered to be composed of a theoretically infinite number 
of small amplitude sine waves of different frequencies. Since 
the wave frequencies differ, the velocities of propagation 
differ (III-l) so that at a particular spot on the ocean the 
phases of the individual waves are uniformly distributed over 
2 7r radians. The wave amplitude at that spot is the sum of 
the individual wave amplitudes and is therefore also random. 

Oceanographers describe the random amplitude of the sea 
by use of an energy density spectrum. This indicates the 
amount of energy in the different component waves which have 
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combined to produce the observed irregular pattern. It is 
a plot of the energy density per unit frequency, S(w) , versus 
frequency, (w) . S(w) is equal to E/ 6 w where E is the energy 
per unit area of the wave surface and 6 w is an incremental 
frequency element. 

The area under the spectrum represents the total energy 
per unit surface area of the sea surface since: 




The area between any two frequencies, w-^, and w^, represent 
the energy density associated with all waves between the two 
frequencies. It can be shown that the total energy in this 
frequency increment 6 w at a central frequency w n is Ref. 

8J7: 

kS(w n ) & w. 

Since the regular waves the energy per unit surface area is 
proportional to the square of the wave amplitude (III-6) it 
follows that the square of the amplitude of a wave having the 
same energy as all the wave components in a band of frequencies 
represented by £ w is: 

A n = 2S(w n } <* w (III-7) 

This allows the construction of an approximate energy density 
spectrum from a finite number of regular wave components for 
a given sea condition. 
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Certain visible properties of the random sea that can be 
measured by trained observers have also been related to the 
total energy E and are listed in Table I /“Ref. 9J7. The 
most commonly referred to property is the Significant Wave 
Height, defined as the average of the heights of the highest 
third of the waves observed and usually written as H-^^* 



Several mathematical models of the energy density spectrum 
have been proposed but the most widely used model today for 
wind generated sea states is the Pierson-Moskowitz spectrum 
/“Ref. 10 _7. It is given by the equation 



S(w) 



. 0.0081 s‘ 



e -0.7Mw o /w)' 



YT 



(III-8) 



. where w Q = (0 . 210 -g/fr^y^) 2 . 

and is plotted in Figure 2 as a function of significant wave 
height. 

Pierson and Moskowitz also developed a relationship 
between the significant wave height and wind velocity, U, 



as : 




0.210 




(III-9) 



TABLE I 

Approximate Wave Height Data 



Most frequent wave height 1.41./E 

Average wave height 1.77 /e 

H i/3 2.83 /e 

Vio 3 ' 60/5 
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Figure 2. Pierson-Moskowitz Energy Density Spectrums 
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Substitution of equation III-9 into III-8 yields 



S(„) = 




e -0.?4gV('»tJ) 4 



(III-10) 



It is this form of the Pierson-Moskowitz spectrum that can 
be used in the XR-3 model to generate wave components for 
the simulation of a random sea as discussed in the follow- 
ing section and Chapter V. 

B. SIMULATION IN THE MODEL 

Two subroutines of the XR-3 Loads and Motions program, 
INCON and WAVES, are involved in the construction and simu- 
lation of waves. INCON develops the wave components to be 
used in the construction of the regular or irregular seas 
and prepares the wave parameter table section of the program's 
summary listing. WAVES calculates the wave forces and moments 
acting on the craft and generates the wave amplitudes at the 
various stations around the hull as well as the bubble volume 
lost due to wave elevation. 

1. Subroutine INCON 

Four options are available to define the wave compon- 
ent frequencies and amplitudes. 

Option 1: The amplitudes and frequencies (rad/sec) 

of up to 10 wave components are specified. 

Option 2: The amplitudes and wavelengths of up to 

10 wave components are specified. 

Option 3 : The significant wave height (H, /,) and 

significant range of periods ar4^specified . 
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Option 4: The significant wave height and range of 

significant wave frequencies are specified. 

Options 1 and 2 allow the user to construct the regular wave 

with one component or irregular seas composed of up to 10 

components. In option 3 and 4 the program calculates the 

wave components from the Pierson-Moskowitz energy density 

spectrum based on the and the range of frequencies in 

the spectrum. If option 3 is used the upper and lower period 

limits specified are converted to the lower (w-^) and upper 

(w u ) frequency limits of the spectrum. 

The desired number of component amplitudes and fre- 
quencies are developed from the spectrum as follows. The 
frequency range is divided into as many increments as there 
are desired wave components. The width of each component is 
w-^c n where n is the increment number and c is a weighted 
scaling factor defined as 

w i/N 
/ u \ 

0 = 

where N is the number of wave components. 

The central frequency of each increment, w n , is deter- 
mined by finding the midpoint of the increment , and the energy 
density per unit frequency, S(w) , is calculated from equation 
III-10 for each w . The central frequencies of each of the 
increments are then used as the component frequencies . The 
component wave amplitudes are found by taking the square root 
of equation III-7, i.e., 

A n = ( 2S(w n ) • 6 w)2 . 
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The above computations are performed in INCON from statement 
INCN3600 to INCN4100. 

After the frequency and amplitude of each wave compo- 
nent is determined, the waveslope, initial encounter frequency, 
and encounter period of each wave component is calculated 
(statements INCN4720 through INCN4-920). This information is 
listed in the wave parameters table section of the program's 
summary listing. It should be stressed that the encounter 
frequencies listed are only true if during the simulation the 
speed is maintained at the initial condition speed. 

2. Subroutine WAVES 

Once the regular wave components with the appropriate 
distribution of amplitude and frequency have been determined 
in INCON, subroutine WAVES constructs the irregular sea simply 
by adding together the regular wave components as time pro- 
gresses. It's important to note that this is not a true 
random sea since the irregular wave formed by the sum of 
several periodic regular waves must also be periodic unless 
of course one of the regular wave frequencies is an irrational 
number. However, it is also apparent that the period of the 
irregular wave may be exceptionally long as compared to that 
of any of the regular waves and that in the course of the 
simulation, one cycle of the irregular wave may not be com- 
pleted. Therefore the irregular sea developed by adding 
together the component regular waves may be considered pseudo- 
random. 
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WAVES generates the wave amplitudes at the various 
stations around the sidewalls and seals relative to calm 
water. The wave elevation relative to any point on the craft 

a. r_ 

is represented for the n x component wave by 

ETA n (x,y) = A sin (-xcosY - y sinY ) + f (t) (ill-ll) 

A n n 

where 

f n (t) = w n t + (-X Q cos/? - Y q si n/3 ) (III-12) 

X Q and Y q represent the position of the center of 
gravity of the craft with respect to the inertial reference 
frame and are time-varying quantities as the craft moves. 

/3 is the angle between the normal to the wave crests and the 
inertial frame X axis as shown in Figure 3* The craft axes 
are x and y and the two define any point on the craft in the 
horizontal plane. The angle Y is the angle between the ship 
x axis and the normal to the wave crests and is defined as 
Y = /3 - V 

It is not difficult to see how the relative motion 
between the craft and wave is accounted for in equation III-ll. 
Assume that the wave height at the center of gravity is desired 
for a craft heading straight into the waves. Then x = 0, 
y= o, /3 = 180° , and H = 0° . Equations III-ll and III-12 
then reduce to one equation 

ETA n (c.g.) = A sin (w t + -^2L X ). 

An u 
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/ 

/ 

wave direction 



Figure J. Relation between ship axes, wave 
direction, and inertial reference 
frame . 
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But X Q is Yt, where V is the craft velocity, and w n = ■ ^ - 

n 

(Eq. III-2) giving 

ETA (c.g.) = A sin( + l-2LX)t. 

A n 7\ n 

The argument of the sine is therefore the encounter frequency 

times time were w = (C + V) as shown in equation A-3. 

e n * 

Equations III-ll and III-12 are represented "by statements 
WAVS2530, WAVS0540, and WAVS0910 . 

Because the integration routine used in the program 
cannot calculate an instantaneous rise in the build-up of 
the wave amplitude, the amplitude of the wave is restrained 
by a factor 

AMPFAC = 1.0 - g-V !• 30 287. 

This restriction of the wave amplitude is evident in the time 
responses shown in chapter 4. 

The primary function of WAVES is to calculate the 
wave forces and moments acting on the craft. This is described 
in detail in reference 1. WAVES also has the option of 
printing the wave elevations at the various stations and 
center of gravity, the plenum volume lost due to the waves , 
and the wave forces and moments. 
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IV. REGULAR SEA RESPONSE CHARACTERISTICS 



A. PROCEDURE 

Determination of the response characteristics of the XR-3 
Loads and Motions computer program was hasically an experimen- 
tal process. In brief, a particular value of craft speed, wave 
amplitude, and wave heading was chosen and held constant while 
the model was repeatedly run in a sinusoidal regular wave, 
each run having a new wave frequency. The desired oscillatory 
response was measured for each run and plotted as a function 
of the encounter frequency. A new value of either craft veloc- 
ity, wave amplitude, or wave direction was chosen and the runs 
were repeated. 

The craft velocities used were 10, 20, and 30 knots which 
covered the speed range of the actual XR-3 test craft. Three 
wave headings of 180, 90, and 0 degrees were chosen thereby 
simulating the craft in ahead, abeam, and following seas 
respectively. The wave frequencies were chosen to cover the 
range of frequencies contained in a random sea and present 
integer values of encounter frequency to the craft as dis- 
cussed in Appendix A. The responses of primary interest 
were roll, pitch, and heave at the center of gravity. 

The XR-3 model can be run with constant speed or constant 
thrust options. With the constant speed option the thrust is 
varied as necessary to maintain the required speed. With the 
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constant thrust option, the thrust is held constant and the 
speed is allowed to vary as the simulation progresses. The 
data presented here was obtained while operating the program 
with the constant thrust option as it simulates the operation 
of the actual XR-3 test craft more closely. However the con- 
stant speed option was used frequently to obtain initial 
operating conditions as is explained later in this chapter. 

The program's initial conditions in pitch, draft, plenum 
pressure, and thrust were determined as a function of speed 
and are given in Appendix C. Appendix D lists the craft 
parameters used. The model was normally run for 35 seconds 
to ensure that the steady-state responses had been obtained. 
The responses were assumed to have reached steady-state when 
their peak-to-peak and average values remained constant for 
at least ten seconds. The responses were normally sampled 
and printed every 0.05 seconds to ensure that the minimum 
Nyquist rate was more than satisfied at the higher encounter 
frequencies. 

1 . Ahead Seas Runs 

The initial runs were performed with a craft speed of 
20 knots in ahead seas ( /3 = 180°) . A wave amplitude of 0.2 
feet was chosen as the excitation and the 20 knot calm water 
thrust of 186.72 lbs. per side was used. Steady-state values 
of heave and pitch were obtained rapidly for encounter fre- 
quencies of 1, 2, and 3 radians per second. However at 
w e = 4 steady-state responses could not be obtained merely 
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by running the model for a greater period of time. As shown 
in Figures 4-A through 4-C, as time progressed the draft 
increased significantly while the pitch oscillations diverged. 
Concurrently the speed had dropped in 60 seconds from an ini- 
tial 20 knots to approximately 15.6 knots thereby reducing 
the frequency of encounter with the wave excitation. More 
important was the fact that the plenum volume and pressure 
were found to have decreased significantly as shown in Figures 
4-D and 4-E. 

Intuitively the above can be explained as follows. As 
the plenum volume and pressure fell, the craft sank deeper 
into the water, increasing the drag on the sidewalls and reduc- 
ing the speed. The only possible source for the loss of pres- 
sure and volume was believed to be associated with the leakage 
of air through gaps formed between the sidewalls and the 
water surface or the bow seal and water surface. It was 
found that the most forward section of the sidewall, and 
therefore presumably the bow seal also, did rise in and out 
of the water allowing air leakage through the gap formed as 
discussed in Appendix E. 

However it should be noted that the pressure-volume 
relationship inside the plenum is an adiabatic process which 
would have the pressure increase instead of decrease for 
constant air mass with a reduced volume. Therefore for 
pressure and volume to both drop, the mass of the air must 
be decreasing at a rate faster than the volume. Since the 
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Figure 4-A. Heave vs. Time 

V=20 kts , A=0 . 2 ft, w =4 

e 

Constant Thrust, (3 =180° 
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Figure 4-B. Pitch vs. Time 
V=20 kts, A=0 . 2 ft, w 0 =4 
Constant Thrust, /3 = 180° • 



36 



25 - 




10 - 



5- 



o i i | ~i 

10 20 30 40 50 

Time (sec) 

Figure 4-C. Speed vs. Time 
V=20 kts, A=0 . 2 ft, w e =4 
Constant Thrust, /2 =180° 
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Figure 4-D. Plenum Pressure vs. Time 
V=20 kts, A=0.2 ft, w q =4 
Constant Thrust, /? =180° 
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Figure 4-E. Plenum Volume vs. Time 
V=20 kts , A=0 .2 ft, w e =4 

Constant Thrust, (3 =180° 
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air flow rate into the plenum will saturate due to the air 
fan map limitations, high leakage rates would cause the mass 
of the air to drop faster than the volume and thereby cause 
the pressure to drop. 

Of interest also is the initial transient portion of 
the responses. As the wave amplitude rises exponentially, 
the responses are almost linear. However as the coupled heave 
and pitch amplitudes increase, gapping begins, resulting in 
the rapid drop in volume and therefore causing the nonlinear 
saturation of the air fans which is reflected in the pressure 
response . 

Clearly the calm water thrust value was not sufficient 
to overcome the increased drag due to the increasing draft. 
Because of this, the model was ran with the constant speed 
option in order to determine the thrust necessary to main- 
tain the 20 knot speed and the encounter frequency of 4 radians 
per second. The result is shown in Figure 5 where the nec- 
essary thrust is a not quite sinusoidal variable with an 
average value of 215 Its. and lagging the wave excitation by 
approximately 240 degrees. The model was then run with a 
constant thrust value of 215 Its. per side and steady-state 
responses were ottained as shown in Figure 6. Note that the 
craft still sinks somewhat deeper into the water than for 
calm water conditions and that the plenum pressure and volume 
have reached a steady-state condition. 
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Figure 5. Starboard Thrust vs. Time 
V=20 kts, A=0 . 2 ft, w g =4 
Constant Speed, /9 =180° 
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Figure 6-A. Heave vs. Time 
V=20 kts, A=0 . 2 ft, w g =4- 
Constant Thrust=215 lb, /3 =180° 
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Figure 6-B. Pitch vs. Time 
V=20 kts , A=0 . 2 ft, w e =4 

Constant Thrust =215 lt>» /S =180° 
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Figure 6-C. Speed vs. Time 
V=20 kts, A=0 . 2 ft, w =4 
Constant Thrust =215 lt>> =180 
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Figure 6-D. Plenum Pressure vs. Time 
V=20 kts , A=0 . 2 ft, w 0 =4 

Constant Thrust =215 lt>» 8 =180° 
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V=20 kts, A=0.2 ft, w g =4 
Constant Thrust =215 lh, /3 =180 



From the initial 20 knot runs a general procedure was 
developed for obtaining steady-state responses. The model was 
initially run with the constant thrust option using the calm 
water thrust for the desired speed. The pitch, heave, and 
roll responses and craft speed were observed. If steady- 
state values were not reached within thirty-five seconds the 
model was rerun with the constant speed option. The average 
value of thrust from this run was then used as the thrust for 
a second constant thrust run. This would normally produce 
steady-state heave and pitch responses for speeds of 20 and 
30 knots. 

Obtaining steady-state responses for a craft speed of 
10 knots required several more iterations of the thrust, how- 
ever. For example, at an encounter frequency of six radians 
per second, the calm water thrust of 200.31 lbs. per side 
drove the craft to a speed of 11. 5 knots in 20 seconds and 
resulted in a decreasing draft and converging heave and pitch 
oscillations. A constant speed run produced an average thrust 
of 19^ lbs. per side but when this was applied as the input 
to a constant thrust run, the speed dropped below hump speed. 

A choice midway between the two values (198 lbs.) again pro- 
duced an increasing speed and unsteady oscillations. A final 
thrust value of 196 lbs. produced the desired steady-state 
responses . 

All runs performed at 10 knots were found to be 
extremely sensitive to small changes in thrust. This is due 
to operating the program in a region where the drag versus 



speed has a minima, allowing small variations in thrust to 
significantly effect the craft speed. This problem was also 
noted by Cagle /~Ref. 11_7 in his study of the original 100-B 
digital model. 

2. Abeam Seas Runs 

To obtain the frequency response functions for the 
roll characteristics of the digital model, the program was 
run simulating the craft in abeam seas, i.e., the angle bet- 
ween the normal to the wave crests (in the direction of pro- 
pagation) and the craft's X Q - axis was 90 degrees. The craft 
was run at 20 and 30 knots using the calm water thrust values 
and the constant thrust option. The wave amplitude was 0.2 
feet. 

A typical set of responses obtained while operating 
the program in this manner is shown in Figures 7-A through 
7-E where neither the heave, pitch, nor roll reached steady- 
state. Figure 7-E provided the explanation for the responses. 
As the craft encountered the waves a yaw oscillation was 
induced. As time progressed the average value of the yaw 
became non- zero and the craft gradually turned into the wave, 
changing the wave excitation-craft response relationship. 

To obtain steady-state pitch, roll, and heave conditions, 
it was clear that it was necessary to first have a stable 
yaw response . 

To obtain steady-state yaw a rudder map, i.e., a set 
of data points that describe the rudder angle as a function 
of time, was added to the model in order to turn the craft 
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Constant Thrust 
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Figure 7-C. Roll vs. Time 
V=20 kts, A=0.2 ft, w 0 =3 
= 90° , Constant Thrust 
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Figure 7-D. Speed vs. Time 
V=20 kts, A=0.2 ft, w g =3 
/3 - 90° , Constant Thrust 
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Constant Thrust 



in a direction opposite that which was induced by the wave. 
Several iterations of this process were required before steady- 
state responses were reached. It was determined that a small 
rudder angle, typically less than a degree, of constant mag- 
nitude would allow steady-state operation to be reached. For 
the example given in Figure 7 , a constant rudder map of -0.2 
degrees as shown in Figure 8-A produced steady-state responses 
as exemplified in Figure 8-B and 8-C. Although the average 
value of the yaw is still not zero, its magnitude is consider- 
ably less than initially when the rudder angle was zero. 

3 . Following Seas Runs 

Following seas are encountered when the craft and waves 
travel in the same direction. The procedures described for 
ahead seas were used to determine the desired response char- 
acteristics . 

B . RESULTS 

The responses of primary interest were heave, pitch, and 
roll. Once steady-state conditions had been reached, the 
peak-to-peak and average values were determined and are listed 
in Appendix G. Since heave is a function of vertical wave 
displacement, the ratio of the two was obtained and plotted 
as a function of encounter frequency. Similiarily the ratio 
of pitch amplitude and roll amplitude to their excitation, 
the wave slope, was found and plotted. 
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!• Frequency Response Functions 



a. Ahead Seas 

Program runs simulating ahead seas conditions 
over a range of encounter frequencies from 1 to 20 radians 
per second were conducted at the following craft speeds and 
wave heights. 

Speed Wave Amplitude 

10 knots 0.2 ft 

20 knots 0.1, 0.2, 0.3, 0.5 ft 

30 knots 0.2, 0.3 ft 



Given that with a sinusoidal excitation a linear system should 
produce a sinusoidal response, in general the heave and pitch 
responses were found to he linear for the ahead seas condi- 
tions. Small nonlinearities were noted at w g =l in the pitch 
response, two examples of which are shown in Figure 9. The 
only extreme occurrence was at V=20 knots, w e =4, and A=0.5 ft. 
As shown in Figure 10 both pitch and heave displayed large 
nonlinear response characteristics. 

Figures 11 through 17 present the frequency 
response functions for heave and pitch for the seven combi- 
nations of craft speed and wave heights simulated. The 
predominate characteristic in all of them is the resonant 
peak in the pitch response near w e =4. A comparison of the 
responses as a function of waveheight for speeds of 20 and 
30 knots is shown in Figures 18 and 19. The heave response 
is nearly independent of waveheight except for the region 
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Figure 9. Ahead Seas Pitch Nonlinearities 
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Pitch Heave 

Figure 10. Ahead Seas Pitch & Heave Nonlinearities 




Figure 11. 



□ = 
0 = 



V=10 kts, A=0 . 2 ft, 
Ahead Seas 

Heave Frequency Response 
Pitch Frequency Response 
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Figure 12 . 

0 = 
0 = 



V=20 kts, A=0 . 1 ft, 
Ahead Seas 

Heave Frequency Response 
Pitch Frequency Response 
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Figure 13. V=20 kts, A=0.2 ft, 

Ahead Seas, 

0 = Heave Frequency Response 
O = Pitch Frequency Response 
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Figure 14. V=20 kts, A=0.3 ft. 

Ahead Seas 

0 = Heave Frequency Response 
O = Pitch Frequency Response 
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Figure 15 . V=20 kts, A=0.5 ft. 

Ahead Seas 

0 = Heave Frequency Response 
G = Pitch Frequency Response 
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Figure 16 . V=30 kts, A=0.2 ft. 

Ahead Seas 

□ = Heave Frequency Response 
O = Pitch Frequency Response 
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Figure 17. V=30 kts, A=0.3 ft, 

Ahead Seas 

0 = Heave Frequency Response 
O = Pitch Frequency Response 
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Figure 18- A. Heave Frequency Response vs. Wave 
Height, V=20 kts 
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Figure 18- B. Pitch Frequency Response vs. Wave 
Height, V=20 kts 
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Figure 19-A. Heave Frequency Response vs. Wave 
Height, V=30 kts 
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Figure 19 -B. Pitch Frequency Response vs. Wave 
Height, V=30 kts 
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near w g =4. This is not true for the pitch response where at 
the higher frequencies the response falls off faster for 
greater wave heights. A comparison of the responses as a 
function of speed for constant waveheight is shown in Fig- 
ure 20. Clearly the response in doth heave and pitch is 
dependent on the speed of the craft. 

A comparison with the four degree-of- freedom 
theoretical response in Figure B-l in Appendix B is also 
of interest. Both models show the speed dependence of heave 
but not of pitch. In general the theoretical response char- 
acteristics fall off slower in both pitch and heave then that 
obtained from the XR-3 digital simulation and the theoretical 
responses would appear to be independent of waveheight at 
resonance. Also noteworthy is the second peak in heave near 
w 0 =8 for the four degree-of-freedom model which does not appear 
in the six degree-of-freedom model until a much higher fre- 
quency is encountered. 

The frequency response functions obtained from 
the XR-3 Loads and Motions program show that at low frequen- 
cies representing waves with large wavelengths as compared 
to the craft length, or the length of the bubble chamber for 
the CAB type craft, the simulated XR-3 will follow the wave 
surface and slope. At high frequencies the wavelength is 
small relative to the plenum length and the craft will tend 
to not follow the waves but will "platform" over them. Inter- 
mediate to these two conditions is a range of resonant fre- 
quencies for which wave inputs, particularly the wave slope, 
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Figure 20-A. Heave Frequency Response Vs. Speed, 
A=0 . 2 ft 
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Figure 20-B. Pitch Frequency Response vs. Speed, 
A-0.2 ft 






7 ^ 



would be amplified to result in craft motions that could be 
significantly larger than the input values. Figure 21 plots 
the dependence of heave and pitch on the wavelength to plenum 
length ratio. Note that in terms of this ratio, heave is 
fairly independent of speed. 

b. Following Seas 

Following seas runs were simulated for craft speeds 

of 10 and 20 knots with a wave amplitude of 0 . 2 feet. Steady- 

state conditions could not be obtained for V=10 knots, w =8. 

e 

Nonlinearities were noted primarily at 10 knots in heave for 
the lower encounter frequencies. Figure 22 presents repre- 
sentative examples. 

The frequency response functions for the following 
seas simulations are presented in Figures 23 and 24. The 10 
knot response characteristics can be compared with the four 
degree-of- freedom theoretical response in Figure B-2. As can 
be seen the two bear little resemblance to each other. In 
particular the magnitudes are considerably smaller for the 
digital simulation. The digital simulation's pitch response 
has the general features of the theoretical response, in 
particular the notch above w g =2 is evident in both, but the 
heave response of the digital simulation has a resonant peak 
in contrast to the notch in the theoretical response. 

c. Abeam Seas 

Runs simulating abeam seas conditions were conducted 
with a wave amplitude of 0.2 feet for craft speeds of 20 and 
30 knots. For abeam seas the frequency of encounter is not 
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Figure 22. Heave Nonlinearities 
Following Seas 
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Figure 23 • 

0 = 
O = 



V=10 kts , A=0 . 2 ft, 

Following Seas 

Heave Frequency Response 

Pitch Frequency Response 
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Figure 24. 



0 = 
O = 



V=20 kts , A=0 . 2 ft 

Following Seas 
Heave Frequency Response 
Pitch Frequency Response 
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a function of the craft speed hut is the actual wave frequency. 
This has the effect of producing considerably greater wave 
slopes than were encountered in the ahead or following seas 
conditions . 

Nonlinearities were found in pitch and roll at 
the lowest encounter frequencies. However since pitch is 
coupled to both the roll and yaw, nonlinearities in its 
response were to be expected. Typical responses are presented 
in Figure 25- 

The heave frequency responses for the two speeds 
were so nearly identical that they can be plotted as one as 
in Figure 26. The frequency response functions for pitch 
and roll are given in Figures 27 and 28 and also show a high 
degree of independence from speed. Comparisons with the 
theoretical responses are not possible since they are based 
on a craft speed of 10 knots. However the theoretical pitch 
response does deserve some attention. The pitch frequency 
response function shows that for the lower encounter frequen- 
cies the pitch amplitude is equal to the wave slope . If 
there were no yaw, as believed to be the case for the four 
degree-of- freedom model, the pitch of the craft would be 
minimal, much as is characterized by the responses from the 
digital simulation. 

2 . Operating Characteristics 

The operating characteristics of the XR-3 program 
were defined as the average draft, average pitch, and the 
thrust necessary to maintain the desired speed. Figures 29 
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Figure 25. Abeam Seas Nonlinearities 




Figure 26. V=20 & 30 kts 

Heave Frequency Response 
Abeam Seas 
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Figure 27. V=20 kts, A=0.2 ft, 

Abeam Seas 

O = Pitch Frequency Response 
A = Roll Frequency Response 



84 





Figure 28. V=30 kts, A=0.2 ft, 

Abeam Seas 

O = Pitch Frequency Response 
A = Roll Frequency Response 
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through 3^ present this information as a function of encounter 
frequency, craft speed, wave height, and craft heading with 
respect to the waves. For both ahead and following seas the 
20 and 30 knot runs show great similarity. At resonance there 
is a large increase in average pitch and draft along with a 
significant increase in the thrust necessary to counter the 
increased drag and maintain the desired speed. However the 
10 knot runs are markedly different in that during resonance 
it was necessary to decrease the thrust to maintain the speed. 

For the abeam seas conditions, the calm water thrust 
was adequate to maintain the desired craft speed. In all 
cases except one, the average draft and pitch remained at 
their calm water values. As noted in the procedures section 
steady-state yaw was a prerequisite for obtaining steady- 
state responses in heave, pitch, and roll. The peak-to-peak 
and average values of yaw at the end of 30 seconds run time 
without rudder control as well as the peak-to-peak value 
with rudder control are listed in Tables G-10 and G-ll of 
Appendix G. 
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Figure 29. Operating Characteristics 
V=10 kts, A=0 . 2 ft, 

Ahead Seas 
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Figure 30. Operating Characteristics 
V=20 kts, A=0.2 ft. 

Ahead Seas 
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Figure 31* Operating Characteristics 
V=20 kts , A=0 . 3 ft, 

Ahead Seas 
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Figure 32. Operating Characteristics 
V=20 kts , A=0 . 5 ft. 

Ahead Seas 
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Figure 33* Operating Characteristics 
V=30 lets, A=0 . 2 ft. 

Ahead Seas 
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Figure 3^* Operating Characteristics 
V=30 kts , A=0 . 3 ft, 

Ahead Seas 
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Figure 35* Operating Characteristics 
V=10 kts , A=0 . 2 ft. 

Following Seas 
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Figure 36 . Operating Characteristics 
V=20 kts , A=0 . 2 ft, 

Following Seas 
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V. IRREGULAR SEA RESPONSE CHARACTERISTICS 



A vessel's behavior in regular waves, as was obtained for 
the digital XR-3 in chapter IV, is of fundamental importance, 
for if the vessel can be described by linear responses, its 
motion in an irregular seaway can be described as the super- 
position of the responses of the vessel to all the wave com- 
ponents of the seaway, i.e., 

Y(w) = X(w) • H(w) 

To check the validity of the ahead sea frequency response 
characteristics obtained for the XR-3 digital model in regu- 
lar seas, the program's ability to simulate the craft's motion 
in irregular seas was used. 

The irregular seas were constructed using wave option 3 
in which given the the ran § e significant periods 

the component wave frequencies and amplitudes are calculated 
by the program using the Pierson-Moskowitz spectrum. Each 
irregular sea generated was composed of ten wave components. 

The significant wave heights and range of significant 
wave periods for various sea conditions were obtained from 
data compiled by the David Taylor Model Basin, a portion of 
which is listed in Table II. The associated spectrums were 
plotted in Figure 2. Using the procedures developed for ahead 
sea runs in regular waves, simulations were conducted for the 
conditions listed in Table III. It was found that in sea 
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TABLE II 



Sea 

State 


Significant 


Wave Heights 


and Wave 


Periods 


Wing 

Velocity 

(knots) 


h/ 3 

(«.) 


%/io 

(ft.) 


Significant 
Range of 
Periods (sec) 


0 


2 


0.08 


0.10 


0-1.2 


1 


5 


0. 29 


0.37 


0.4-2. 8 


1 


8.5 


1.0 


1.2 


o .8-5.0 


2 


10 


1.4 


1.8 


1. 0-6.0 


2 


12 


2.2 


2.8 


1. 0-7-0 


2 


13.5 


2.9 


3-7 


1.4-7. 6 



Irregular Sea 


TABLE III 

Simulation Conditions 


Thrust 
per side 
(16) 


Craft 

Speed 

(knots) 


R l/3 

(ft.) 


20 


0.29 


186.72 




1.0 


210.0 




1.4 


223.0 


30 


1.0 


300.0 




1.4 


310.0 
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states with significant wave heights greater than those listed 
in Table III there was nearly continuous contact with the top 
of the plenum chamber, a condition which is not accounted for 
in the program and therefore avoided. A typical example of 
the generated irregular wave as seen at the craft's center of 
gravity and the responses to the wave are shown in Figure 37-A 
through 37-E. Clearly apparent in the heave and pitch responses 
is the "low pass filtering" of the wave excitation's high 
frequency components which was the dominate characteristic 
of the heave and pitch regular sea frequency response functions. 

The magnitude of the frequency response function is 



H(w) 




In order to determine the excitation and response spectrums, 
the discrete fourier transform (DFT) was taken of the wave 
excitation using an IBM 3^0 library subroutine as described 
in Appendix F. 

To obtain an accurate representation of the spectrums 
with the DFT, several conditions had to be met. To prevent 
aliasing, the high-frequency components of the time function 
impersonating low frequency components, the sampling rate of 
the time response or excitation had to be greater than the 
Nyquist rate. A sampling interval of 0.03 seconds satisfied 
this condition. The frequency resolution of the DFT is l/c 
where is the sample length. The longer the sample the 
greater the resolution. However this conflicted with the 
actual running of the program in which greater than 2.5 minutes 
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Wave Height (ft) 



0 , 6 - 




- 0.6 - 



Figure 37-A. C.G. Wave Height vs. Time 
V = 20 kts , H l/3 = 1.0 ft 
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Heave (in) 




Figure 37-B. Heave vs. Time 
V = 20 kts, ^ 1/3 = 1«0 ft 
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Pitch (degrees) 




Figure 3 7 -C. Pitch vs. Time 
V = 20 kts , H l/3 = 1.0 ft 
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300 A 
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20 kts 



Pressure (lb) 




Figure 37-E. Plenum Pressure vs. Time 
V = 20 kts, H 1 = 1.0 ft 
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CPU time was required to simulate the craft for one second of 
run time in irregular seas. Finally, the bandwidth of the 
spectrum is a function of the number of sample points and 
the sample length /“Ref. 12 \_J , A 512 point DFT, which there- 
fore produced 256 frequency components, and a sample length 
of approximately 1 6 seconds were found to produce good discrete 
spectrums . 

Figure 38-A shows the discrete magnitude spectrum of the 
wave height at the center of gravity for V=20 kts, 1^=1. Oft. 
Note that the spectrum produced is a function of encounter 
frequency since the data was taken while the craft was in 
motion. Also shown in the figure is the encounter frequencies 
and amplitudes of the components used in actually constructing 
the irregular wave. The two agree quite well. 

Figures 38-B and 38 -C present the magnitude spectrum of 

the heave and pitch responses to the above wave excitation, 

again as a function of encounter frequency. The "low pass" 

characteristics of the heave and pitch regular sea frequency 

response functions are clearly evident in the spectrums also. 

However also noted are the components at frequencies below 

w =3 for which there is no excitation. These values were 
e 

ignored. From the wave spectra data, the frequencies and 
amplitudes of the peaks were obtained. The amplitudes of the 
heave spectra were also obtained at the same frequencies and 
the ratio of the two were calculated (see Table G-12, G-13* 
and g- 14) and are plotted in Figure 38 -D. As a comparison 
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Wave Amplitude (in) 




Figure 38 -A. C.G. Wave Height Spectrum 
V= 20 kts, = 1.0 ft 
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Heave Amplitude (in) 




Figure 38 -B. Heave Spectrum 

V = 20 kts, H 1 = 
105 



1.0 ft 



Pitch Amplitude (degrees) 




Figure 38-C. Pitch Spectrum 

V = 20 kts, = 1.0 ft 
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Figure JQ-B. Heave Frequency Response 
V=20 kts, 3-^=1. 0 ft 



107 




Figure 38 -E. Pitch Frequency Response 
V=20 kts, R^yl.O ft 
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with the heave frequency response resulting from regular seas, 
the response for V=20, A=0.3 ft. is also plotted. 

To determine the irregular sea pitch frequency response 
it was necessary to first calculate the actual wave frequency 
from the encounter frequency and then the maximum wave slope 
for each wave frequency. The ratio of the pitch amplitude to 
the maximum wave slope was then determined and is plotted in 
Figure 38-E. Note that for both the heave and pitch irregular 
sea frequency response functions there is close agreement 
with those obtained in the regular seas runs. 

The spectrums of wave amplitude and pitch and heave 
responses for two other combinations of sea state and craft 
speed are shown in Figures 39 through 40 . Also presented 
are the associated frequency response functions. As can be 
seen, the frequency response functions for the irregular sea 
runs agree well with the regular sea responses. 
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Wave Amplitude (in) 



3.0 - 



2.5 - 



2.0 - 



1,5 - 



l.o - 




Figure 39-A. C.G. Wave Height Spectrum 
V = 20 kts , = 1.4 ft 
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Heave Amplitude (in) 




Figure 39-B. Heave Spectrum 

y = 20 lets, h i/3 = 1.4 ft 
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Pitch Amplitude (degrees) 




Figure 39-C. Pitch Spectrum 

V = 20 kts , = 1.4 ft 
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Figure 39-D. Heave Frequency Response 
V=20 kts , . k ft 
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w. 



Figure 39-E. Pitch Frequency Response 
V=20 kts , 
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Wave Amplitude (in) 




Figure 40-A. C.G. Wave Height Spectrum 
V = 30 kts , Rwj = 1.0 ft 
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Heave Amplitude (in) 




Figure 40-B. Heave Spectrum 

V = 30 kts, H 1 = 1.0 ft 
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Pitch Amplitude (degrees) 



0.7 




Figure 40-C . Pitch Spectrum 

V = 30 kts, = 1.0 ft 
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w. 



Figure 40-D. Heave Frequency Response 
V=30 kts , H 1//3 =1.0 ft 
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w. 



Figure 40-E. Pitch Frequency Response 
V=30 kts, K^yl.O ft 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The XR-3 Loads and Motions Program was subjected to a 
wide variety of wave conditions. The conditions necessary 
to obtain steady-state operation were determined and the 
frequency response functions for heave, pitch, and roll were 
developed. The heave and pitch frequency response functions 
obtained from the regular ahead seas runs were found to com- 
pare well with those obtained from irregular ahead seas by 
use of the DFT. In general the nonlinearities inherent in 
a six degree-of- freedom model were found to be minimal and 
exist at the lower encounter frequencies. It is believed 
that the results from this study provide a good baseline 
with which to compare future modifications to the model. 

A most interesting study would be to obtain frequency 
response data from the actual XR-3 test craft and compare 
it with that obtained from the digital model. It is recom- 
mended that a closer look be given to the model's simulation 
of the craft at the lower speeds. It was noted that during 
resonance at 10 knots, as the craft sank deeper actually 
less thrust was required to maintain the desired speed than 
at the other encounter frequencies. There appears to be no 
plausible explanation for this. 

Although the conditions under which gapping of the side- 
walls occurred were determined no effort was made to determine 
the actual air flow rates out of the gaps or if indeed the 
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air loss was great and rapid enough to cause the loss in plenum 
pressure noted and therefore the settling of the craft deeper 
into the water. Finally, although not analyzed or presented 
in this study, data was obtained and is available on the 
vertical accelerations of the center of gravity for the major- 
ity of the steady-state runs. Analysis of this data could be 
of use in defining the habitability zones for the simulated 
XR-3. 
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APPENDIX A 



Encounter Frequency 



The frequency of excitation is commonly given in terms 
of the encounter frequency, the frequency actually "seen" by 
the vessel. Due to the forward speed of the craft, the fre- 
quency of encounter with the waves is increased or decreased 
depending on the relative direction between the vessel's 
course and the direction of propagation of the surface wave. 
The encounter frequency is given by: 



w g = ~^ L - (C-V cos (i ) (rad/sec) 



(A-l) 



where = wave length (ft) 



C = wave velocity (ft/sec) 

V = craft velocity (ft/sec) 

/3 = the angle between the craft’s course and 



the direction of propagation of the wave. 



Given that the wave velocity, C — ( 5 — — ) 2 > and the wavelength, 

^ IX 



where w is the wave frequency, by substitution into 



w 

A-l it can be shown that: 



w = w-w • 
e 



2 V cos 3 . 
g 



(A- 2) 



For a vessel heading into the waves, /3 = 180° and 



(C+V) . 



(A-3) 
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For this study it was desired to simulate regular waves 
that produced integer values of w g between 1 rad/se'c and 20 
rad/sec for a given value of vessel velocity. The application 
of the quadratic formula to equation A-2 yields 



w 



1 - J 1+ (4Vw e )/g 



-2 




(A-4) 



for (3 = 180°. 

This relationship is plotted in Figure A-l. 
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Figure A-l. Wave Frequency vs. Encounter 
Frequency 
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APPENDIX B 



Theoretical Frequency Response Functions 

A theoretical frequency response model for the XR-3 was 
developed by the Aerojet-General Corp. /~Ref. 7_7 using a four 
degree-of-freedom representation of the craft with responses 
in heave, pitch, roll, and surge. It was assumed that the 
response motions were small enough so that deviation from 
the linear response was negligible for the motions in the 
vertical plane. The requency response characteristics obtained 
are reproduced in Figures B-l thru B-3 and show not only the 
theoretical response but the results from |r scale model tow- 
ing tests as well. 
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Figure B-1A. Heave Frequency Response, 
Ahead Seas 




Figure B-1B. Pitch Frequency Response, 
Ahead Seas 
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Figure B-2A. Heave Frequency Response, 
Following Seas 




U> e ^ TMCCOMTCC. FtecQuCMOr (:zo" 

Figure B-2B, Pitch Frequency Response, 
Following Seas 
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Figure B-3A. Heave Frequency Response, 
Oblique Seas 




Figure B-3B. Pitch Frequency Response, 
Oblique Seas 
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Figure B-JC. Roll Frequency Response, 
Oblique Seas 
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APPENDIX C 



Initial Conditions 

Prior to running the computer model with simulated wave 
conditions the "calm water" initial conditions had to be 
determined. Previous work by Forbes /~Ref. 4_J7 and Menzel 
/~Ref. 5_7 had shown that the initial conditions used in the 
program had to be close to the steady state operating values. 
Failure to ensure this would result in either an imbalance in 
the forces generated within the program and cause the program 
to stop due to an improper integration step size or transients 
that obscured the response of the variable under study. 

The steady state conditions found by Menzel for his 
Program 3 /Table XVI, Ref, 5_7 were used as the initial input 
to the computer model. The program was run using the constant 
speed option for forty seconds simulation time at nine different 
speeds. The steady state value of pitch angle, draft, plenum 
pressure, and thrust were determined for each speed. These 
values were then used as the initial conditions for a second 
set of constant speed runs to verify that the steady state 
values produced no transients when used as initial conditions. 

A third set of runs with the constant thrust option in use 
was performed to verify that the speed was maintained at its 
desired value. 
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The "calm water" steady state values determined by the 
above procedure and used throughout this study as initial 
conditions are listed in Table C-l and plotted in Figure C-l. 



TABLE C-l 





Calm Water 


Initial 


Conditions 




Speed 


Pitch 


Draft 


Plenum 


Thrust 




Angle 




Pressure 


per sid< 


(knots) 


(deg) 


(in) 


(psf) 


(lb) 


10.0 


1.62 


8.17 


23-93 


200 .31 


12.5 


1.11 


7.03 


24.84 


174.71 


15.0 


0.84 


6.74 


24.84 


167.85 


17.5 


O .63 


6.41 


24.84 


173.23 


20.0 


0.48 


6.12 


24.84 


186.72 


22.5 


O .36 


5-87 


24.84 


205.77 


25.0 


0 .29 


5.66 


24.84 


229.31 


27.5 


0.25 


5.48 


24.84 


256 . 65 


30.0 


0.26 


5.34 


24.84 


287.22 



131 




132 




APPENDIX D 



Craft Parameters 

The following craft parameters were used throughout this 
study. The parameters are listed by Input Blocks. Also 
listed are the associated variable names as used in the 
program. 



Integrator Tolerances: Block 1 

Number of Integrators Used: 10 

Runga-Kutta Variable Step Sized Used 

-4 -4 -5 -7 

Initial Step Sizes: 1x10 , 2x10 , 1x10 , 1x10 , 

lxlO" 4 , lxlO -6 , lxl0~ 9 , lxlO -9 , lxlO'\ lxlO -4 



Masses and Inertias: Block 2 

Weight: 6050 lb WEIGHT 

Longitudinal Center of Gravity: 10.05 ft forward of 

transom, XS 

Vertical Center of Gravity: 2.5^ ft above keel, ZS 

2 

I , Mass Moment of Inertia About X axis: 2870.0 slug-ft 

xx 



AIXX 

I , Mass Moment of Inertia About Y axis: 

yy 

AIYY 



9320.0 slug-ft 2 , 



I Mass Moment of Inertia About Z axis: 10,580.0 slug-ft 

zz 

AIZZ 

I , Mass Moment of Inertia About XZ axis: -2800.0 slug-ft 

xz 



AIXZ 



Craft Geometry: Block 3 

Total Craft Length: 24. 7 ft, XLTOT 

Number of Stations along Port and Starboard Sidewalls: 11 
Number of Stations along Bow and Stern Seals : 5 

Sidewal: Block k 

Y distance from centerline to sidewall: 5*37 ft, YSW 

Average Wetted Length of Sidewall: 21.9 ft, XLSW 

Leakage Orifice Coefficient of Sidewall: 0.7, CFSW 

Cross-flow Drag Coefficient of Sidewall: 1.28, CDSW 

Average Beam of Sidewall: 0.5 ft, AVBMSW 

Stern Seal: Block 5 

X corrdinate of Stern Seal Hinge: 3*97 ft forward of 

transom, XSSI 

Z coordinate of Stern Seal Hinge: 1.875 ft above keel, 

ZSSI 

Base leakage area: 3*87 ft^, ALEAK 

Seal Leakage Orifice Coefficient: 0.9. CFSS 

Length of Rear Support Cable: 1.875 ft 

Pressure Differential Between Stern Seal and Plenum 
Chamber: 1 psf, DPSS 

Length of Leading Edge of Seal: 4.06 ft, XLF 

Bow Seal: Block 6 

X coordinate of Bow Seal Hinge: 23-44 ft forward of 

transom, XBSI 

Seal Leakage Orifice Coefficient: 0.9, CFBS 

Pressure Differential between Bow Seal and Plenum 
Chamber : 1 psf , DPBS 
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2 coordinate of Bow Seal Hinge: 1.875 ft above keel, 2BSI 

Length of Rear Support Cable: 1.875 ft 

Length of Leading Edge of Seal: 4.0 6 ft, XBF 

Base Leakage Area: 0.08 ft 2 , BLEAK 

Length of Middle Support Cable: 1.875 ft 

Plenum: Block 7 

Plenum Length at Water Surface: 20.0 ft, XLBW 

Plenum Width at Water Surface: 10.0 ft, XBBW 

X coordinate of Pressure Wave Pivot Point: 17*2 ft forward 

of transom, XFMV 

Plenum Length at Deck: 20.0 ft, XL 

X coordinate of Center of Pressure: 10.4 ft forward of 

transom, XCP0 

Plenum Average Height: 1.915 ft, BUBHGT 

Froude Number: 0.556, FNCRIT 

Propulsion Geometry: Block 8 

X coordinate of Propeller Center: -1.275 ft forward of 

transom, XPO 

Y distance from centerline to propeller center: 5-55 ft, 

YPO 

Z coordinate of propeller: -0.604 ft above keel, ZP0 

Rudder Geometry: Block 9 

X coordinate of Centroid of Rudder, -1.125 ft forward 
of transom, XR0 

Y distance from Centerline to Rudder Centroid: 5-55 ft. 

YR 
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Z coordinate of Centroid of Rudders -0.208 ft above keel, 
ZRO 

Rudder Spans 1.21 ft, RSPAN 
Rudder Aspect Ratios 2.15, RASPR 
Rudder Areas 0.68 ft 2 , RAREA 

Average Thickness Ratio of Rudder Sections 0.167, RTC 
Aerodynamics s Block 10 

Reference Lengths 20.0 ft, XLAER0 
Reference Widths 10.0 ft, BEAM 
Thrust: Block 16 

Starboard and Port Side Thrust Coefficients s 0.01, STHS 
and STHP 
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APPENDIX E 



Gapping 

Subroutine SIDEWL calculates the forces and moments acting 
on the craft due to the sidewalls as well as the leakage flow 
rates associated with any gaps which open under the sidewalls 
due either to craft motion or waves. A print switch option 
allows printing the sidewall gaps, imersion depths, and total 
forces and moments at each integration step. For purposes 
of this study in which 35 second runs were normally conducted 
the printed output of this print option was excessive since 
all the above quantities are printed each integration step. 
Since it was desired to know only if a gap occurred, its 
location along the sidewall, and the magnitude, a simple 
program addition was made to SIDEWL as follows. 

C GAP OR WETTED DRAFT CALCULATION SDWL0460 

C SDWL0^70 

*- GAPT=0 . 0 

DO 10 J=1 , 2 SDWL0480 



7 GAP(J,K)=-DDIN*(1.-(DSW(J,K) )/PBHEAD 
GAPT=GAPT+GAP( J,K) 



GO TO 10 
8 GAP( J,K)=0 .0 
10 CONTINUE 



< 



IF( GAFT .EQ . 0 . 0 ) GO TO 12 
WRITE (6,U)VAL(1) ,((GAP(I,J),J=1,H) » 
11=1,1) , GAPT, ((GAP(I.J), 1=1,11,1=2,2) 
11 F0RMAT( ' (T) ' ,F10 .5,11(F8.5) ,/, 

2’(G)' ,FL0.5*H(F8.5)) 



SDWL0630 



SDWL0640 

SDWL0650 

SDWLO 66 O 
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SDWL 0670 

SDWL0680 



C 

C LEAKAGE AREA 

—*"12 ALSW= 0 . 0 

Gapping was found to occur under the sidewall station 1 
for the following conditions: 

V = 10 knots, A = 0.2 ft 

w = 4 
e 

V = 20 knots, A = 0.2 ft 

w q =4 . 5 

A = 0.3 ft 

w 0 =4, 5 , 6 , 8 , 20 

A = 0.5 ft 

w 0 = 4, 5 , 6 , 8 , 10 

V = 30 knots, A = 0.2 ft 

w 0 = 4, 5, 6 , 8 , 10 

A = 0.3 ft 

w e = 4, 5, 6 , 8 , 10 , 15 

Gapping was found under station 2 for the following conditions: 

V = 20 knots, A = 0.5 ft 

w g = 5, 6, 8, 10 

V = 30 knots, A = 0.3 ft 

w 0 = 5 , 6 , 10 , 15 

Gapping was noted at additional stations for V=20 knots, 

A = 0.5 ft. No gapping was noted during steady-state operation 
in abeam seas, however it was noted under the most forward 
station in following sea runs at V=10 , w g =4 and at station 
2 at V=20, w 0 =4. 
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APPENDIX F 



A DFT Program 

The IBM 360 library subroutine HARM which performs 
discrete complex fourier transforms on a complex three 
dimensional array, was used to take the DFT of the wave 
excitation and the heave and pitch time responses from 
ahead irregular seas in order to determine the frequency 
response functions. Once the complex fourier transform 
of the excitation and responses was determined, the 
magnitude of the spectrum was computed and plotted using 
the library subroutine DRAWP. The program is listed on 
the following pages . 
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c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 



c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 



C I M ENS I CN W ( 600 ) , Z ( 600 ) » P ( 600 I 

CIMENSICN FV ( 514 l 

CIMENSICN S(1012),INV( 1012) »M ( 3 ) 

CIMENSICN TIM E( 600 ) » WH ( 600 ) , Z0 ( 600 ) , PA ( 6C0 ) 
CCMFLSX*8 A ( 1024 « 1 , 1 ) 

INT EGER*4 IT8 ( 1 2 J / 12*0/ 

REAL*4 PT8<28 )/26*0.0/ 



REAC IN THE POWER CF TWC (DETERMINES THE NUMBER 
CF CATA POINTS ) SAMPLE LENGTH, ANC SAMPLING 
INTERVAL 



REAC(5,5)KK.T,DT 

FCPMAT( I5,F10.4,F10.4) 

P 1=2.14 15S 

N1=2**KK 

N 2 s 1 

N 2 = 1 

NF=Nl/2 



REAC IN TIME, WAVE HEIGHT, HEAVE, ANO PITCH CATA 



CC 10 K = 1 »N1 

REAC (5,6 >TIME(K) ,WH(K I ,ZD(K ) ,PA (K) 
6 FCPMAT ( 4F16.4) 

10 CONTINUE 
M <1) = KK 
M ( 2)=0 
M (2 )=0 



CALCULATE THE FUNCAMENTAL FRECUENCY ANO THE 
FPECUENCY COMPONENTS 



FLNC=(2.0*PI )/T 
F V ( 1 )=0 *0 
CC 20 J = 2 , NH 
FV( J)=FV (U-l )+FUND 
20 CONTINUE 



SET THE COMPLEX A MATRIX TO ZERO 



CC 2 11=1, N1 
CC 2 12=1, N2 
CC 2 13=1, N3 
2 A(I1, 12, I3)=(0. 0,0.0) 



TRANSFER THE WAVEHEIGHT DATA TO THE COMPLEX A MATRIX 



CC 2 11=1, N1 
CC 3 12=1, N2 
CC 2 13 = 1, N3 
A (11,1, 1 )=WH(I1) 

2 CONTINUE 

COMPUTE THE DISCRETE FCLRIER TRANSFORM CF THE 
V A V E HEIGHT DATA 

CALL HARM(A,M, I N V ,S ,- 1 , 1 FSRR ) 
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C 

C 



c 

c 

c 

c 



c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 



c 

c 

c 

c 



CALCULATE THE MAGNITUCE OF THE FOURIER COMPONENTS 
ANC CONVERT THE RESULTS TO INCHES 



CC 20 11=1, NH 

20 W (I1)=CABS(A(I1, 1, 1) )*24.G 
CLEAR THE A MATRIX 



CC A 11=1, N1 
CC 4 12=1, N2 
CC 4 12=1, N3 
4 A (11,12 , 13 ) = < 0. 0,0.0) 

COMPUTE THE D FT OF THE HEAVE CATA AND THE MAGNITUDES 
OF THE FOURIER COMPONENTS 



CC 7 11=1, N1 
CC 7 12 = 1, N2 
CC 7 12 = 1, N3 
A (II, 1, 1 )=ZD(I1 ) 

7 CONTINUE 

CALL HAPM(A,M, I NV , S ,-l , I FERR ) 
CC 50 1 = 1, NH 

50 Z(I)=2.0*CABS(A(I,1,1) ) 



THE CC VALUE OF HEAVE (THE AVERAGE CRAFT) IS 
STORED AND LATER FRINTEC. IT IS NCT PLOTTEC 



ZAVG=Z( 1 ) 
Z(1)=0.C 



CCMFUTE THE DFT OF THE PITCH ANGLE CATA ANC THE 
MAGNITUDES OF THE FOUR I ER COMPONENTS 



CC £ 11 = 1, N1 
CC 6 12=1, N2 
CC e 13=1, N3 

8 A(I1,I2,I3)=(0.0,0.0) 

CC 9 11=1, N1 

CC 9 12 = 1, N2 
CC 9 12=1, N3 
A(I1, 1,1)= PA (II) 

9 CONTINUE 

CALL HAFM(A,M,INV,S,-1,IFERR) 
CC 55 J = 1 , NH 

55 F U)=2.C*CABS( A( J,l,l) ) 



STORE THE CC VALUE CF PITCH FOR LATER OUTPUT. IT 
IS NOT FLOTTED 



F AVG = P( 1 ) 

F (1 ) = 0.0 

W FITE (6,60 ) T,DT, M,F UNO, ZAVG, PA VC- 
tO FCFM AT ( / , ' SAMPLE TIME = ',F7.4,' SECONDS'//,' SAMP 
SLE INTERVAL = *,F6.4,' SECGNCS'//,' NUMBER CF SAMFL 
€ E S = ',14//,' FUNDAMENTAL FREQUENCY = *,F6.4,' RPS' 
£//,' ZAVG = ' , F7 • 4// , * PA VG = *,F7.4) 



PLOT THE SPECTRUMS 
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I T E { 12 ) =1 
I7E (3 )=5 
PTE (1)=2Q*C 
PTE (2)=0.5 

C/ll CRAIrtP (NH, F V » b » I TB »RTB ) 
PTE (1 )= 10.0 
P TE ( 2 J = 0 .2 

C/ll CP/WP (NH, FV » Z » I TE »RTB ) 
PTE (1 )=E.O 
RTE (2) = C.l 

C/LL CR/WP (NH,FV,PtITB t RTB) 

STCP 

ENC 
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APPENDIX G 



Data Tables 

The frequency response and operating characteristics data 
presented graphically in the text are listed in Tables G-l 



through G-14. 


The undefined headings and symbols are : 


Ratio .1 


the ratio of the peak-to-peak heave to the 
peak-to-peak wave height 


Ratio 2 


the ratio of the peak-to-peak pitch amplitude 
to twice the maximum wave slope 


Ratio 3 


the ratio of the peak-to-peak roll amplitude 
to twice the maximum wave slope 


A 


wavelength (ft) 


L 


plenum chamber length, 20 ft. 


* 


water contact with the top of the plenum 
chamber occurred during steady-state operation 


Z 


calm water draft (in) 


9 


calm water pitch (degrees) 


T 


calm water thrust used per side (lb) 


V 


craft velocity (knots) 


A 


wave amplitude (ft) 


0 


angle between the craft and wave headings 
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